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Abstract  Space-variant  ultrasound  attenuation  in  tissue  is  one  of 
the  phenomena  that  should  be  taken  into  account  when 
reconstructing  images  from  ultrasonographic  data.  This 
contribution  presents  a  method  to  estimate  automatically  one¬ 
dimensional  attenuation  profiles  along  isolated  “rays”,  forming 
ultrasonic  images.  The  estimates  are  based  entirely  on  digital 
ultrasonographic  data,  without  any  reference  to  a-priori 
knowledge  on  anatomical  structures  involved.  Using  a  standard 
model  of  ultrasound  propagation,  the  “layered”  and  discretised 
formulation  allowed  for  simplified,  thus  manageable,  least- 
squares-error  estimation  of  the  attenuation  profiles.  The  results 
of  analysing  ultrasonographic  data  obtained  from  biological 
objects,  as  presented  on  figures,  are  rather  promising. 
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I.  Introduction 

Contemporary  medical  B-mode  imaging  uses  wide-band 
ultrasonic  impulses  that  are  attenuated  due  to  several  physical 
phenomena,  when  travelling  through  the  imaged  tissues.  The 
space-variant  and  frequency-dependent  attenuation  should  be 
explicitly  taken  into  account  when  reconstructing  the  image 
to  prevent  disturbing  artefacts.  Several  papers,  e.g.  [3-5], 
have  been  published  on  methods  of  reconstructing  the  image 
using  a  known  attenuation  map  of  the  imaged  area,  but  a 
reliable  and  simple  estimation  of  the  attenuation  distribution 
is  still  a  challenging  problem.  When  excluding  the 
approaches  based  on  a-priori  known  tissue  types  in  the  range 
of  the  image,  the  attenuation  distribution  has  to  be  estimated 
automatically  by  analysing  the  received  standard  ultrasonic 
echo,  as  any  other  measurements  are  usually  not  feasible  in 
clinical  routine.  The  presented  contribution  introduces  a 
simplifying  approach  enabling  to  determine  individual 
attenuation  profiles  along  each  image  ray  from  the  measured 
ultrasonographic  data  automatically,  taking  into  account  the 
properties  of  ultrasound  propagation. 
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II.  Methodology 

The  sonographic  scan-head  transmits  a  wide-band 
ultrasonic  impulse  that  travels  along  a  straight  path  (the 
“ray”)  with  the  decreasing  wave-amplitude  as 

Azf)=A0,f)e^J’z)z,  (1) 

where  A(z,  f)  is  the  amplitude  of  the  wave  component  at  the 
frequency  f  in  the  distance  z  from  the  transducer  and 
ju(f,  is  the  local  attenuation  coefficient.  Let  us  consider  a 
certain  frequency  range  of  the  signal,  f  e  (  f  ,  /2 ) ,  and  the 
propagation  with  only  a  single  reflection  or  scattering  on  each 
inhomogeneity  on  the  path,  thus  neglecting  any  secondary 
reflections.  The  ultrasound  absorption  coefficient  is,  as  in  [6], 
supposed  linearly  dependent  on  frequency,  ju(f)=  p  f  ■ 
Constant  velocity  c  of  ultrasound  propagation  in  tissue  is 
assumed. 

The  used  model  of  the  imaged  tissue  is  formed  as  a 
composition  of  M  homogeneous  layers  along  a  ray  (fig.  1). 
Every  layer  has  its  constant  attenuation  coefficient  p. ,  its 

depth  A,  and  its  local  signal-retum-ratio  o  ,  produced  both 

by  specular  reflection  and/or  backscattering  distributed  in  the 
layer.  If  there  are  specular  reflections  due  to  abrupt  acoustical 
impedance  changes  on  the  ray  trace,  the  corresponding  thin 
areas  are  treated  as  separate  layers. 
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The  transducer  emits  a  short  impulse  at  time  t  =  0  and 
then  it  scans  the  returning  echo.  The  signal  received  in  time 
/  was  obviously  reflected/scattered  in  the  tissue  in  time  '2 1] 

and  it  has  passed  M(/,)  layers,  each  of  them  twice,  the 

u  (/,  )-i 

M  (/j)  -th  layer  only  to  the  depth  ±^c  _  ^  A .  •  ^  can  easily 

/= i 

be  shown,  based  on  (1),  that  the  amplitude  of  the  echo 
component  at  the  frequency  f  is 


A.K,f)  =  Ad  /2,/)  =  ^0,/)exp  -/2]jUA-  Gm(,Y  (la^ 


«(>) 


On  the  first  look,  it  may  be  argued  that  this  expression  does 
not  take  into  account  the  phase  relations  among  individual 
returned  components  of  the  wave.  Nevertheless,  according  to 
the  problem  definition,  we  are  not  trying  to  compensate  for 
the  interference,  causing  e.g.  the  speckles,  but  exclusively  for 
the  differing  attenuation  of  the  signal  received  at  different 
time  instants. 


Applying  logarithmic  transform  to  (la),  which  is  always 
possible  as  all  its  components  are  positive,  yields 

ln^/)  =  ln^0>/)-/X>uA  +%(,) 

i= 1 

for  the  interval  (/  ,f2).  Here,  s  can  be  interpreted  as  a 

logarithmic  measure  of  the  local  signal-retum-ratio.  The 
primary  aim  of  the  analysis  is  to  determine  the  vector  jj.  but 
as  the  only  quantity  accessible  to  measurement  is  A&f ),  all 
other  quantities  on  the  right  hand  side  of  (2)  must  be 
determined  or  estimated  as  well. 


When  using  the  (Hilbert-transform  derived) 
envelope  of  the  received  wide-band  radio-frequency  signal, 
we  obtain  the  resulting  integral  amplitude  of  the  echo,  which 
can  be  regarded,  on  average,  approximately  proportional  to 
the  amplitude  mean  value  with  respect  to  frequency.  Working 
with  the  mean  values  removes  the  frequency  dependence 
from  the  formulation.  Although  this  approximation  may  seem 
to  ignore  one  important  aspect,  the  experimental  results 
support  it  as  a  sound  simplification,  making  the  problem 
reasonably  tractable.  This  way,  we  obtain 


M(t) 


+  s  > 


/A 


w- 


fi 

Let’ s  denote  y[z)=  jin  A(z,f)df  and  approximate  it  by  the 

A 

measured  signal  value  .  After  a  few  simplifying  steps 

and  integration  we  get 


Our  aim  is  to  estimate  Ju/  and  s .  numerically,  utilising 
the  known  measurement  j(^).  The  tissue  model  is  initially 
chosen  as  a  (small)  number  of  layers  with  corresponding 
thicknesses  A.  and  is  then  automatically  refined  during  the 
computation  thus  enabling  for  arbitrarily  complex  structure. 
For  the  computation,  the  equation  (3),  after  some 
modifications,  is  discretised  in  z  with  the  space 
sampling  period  A  %  =  cT  ,  where  the  time  sampling  period 
T  is  given  by  the  A/D  converter  of  the  scanner.  The  total 
number  N  of  samples  on  the  ray  is  therefore  fixed  while  the 
number  M  of  layers  approximating  the  tissue  structure  is 
iteratively  updated.  After  discretisation,  the  mathematical 
model  acquires  the  matrix  form 

Xu  =  1 ,  (4) 

where 


1  = 


M) 

J 1  -Jo  +/3  2-XA  f 4^  0,U  (l) 
;= l 

M  (N)  ~ 
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(u  containing  the  sought  corrections  d/J  and  dsr  to 
arbitrarily  chosen  initial  estimates  /uoi  and  j-  see  [9]) ,  and 
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Because  the  size  of  X  is  Ax 2 M  ( N >  2 M  ),  the 
equation  system  is  over-determined  and  can  only  be  solved  in 
the  sense  of  square-error  minimisation.  According  to  [8],  this 
estimate  is 


TT  =  (x7x)~‘xr  I.  (4a) 


M(,) 


j(^)=j(0)-y(/22  +(/2  ~fi)s 


(3) 


The  details  of  the  iterative  algorithm  including  the  initial 
estimates  and  the  refinement  strategy  are  described  in  [9], 


III.  Experimental  results 

The  designed  algorithm  has  been  extensively  tested  on 
radio-frequency  ultrasonographic  data  measured  on  both,  a 
well  identifiable  tissue-phantom,  and  a  biological  in-vivo 
specimen  (pig’s  heart).  The  data  were  scanned  with  2964 


Fig.  2  Measured  phantom  A-mode  signal  (detected  envelope) 


radio  frequency  samples  on  a  ray,  in  the  frequency  range 
from  2.5  to  7.5  MHz,  and  the  number  of  rays  in  one  frame  (a 
fan)  was  78. 

Example  results  of  the  extensive  experimentation  are 
presented  on  figures.  On  fig. 2,  the  Hilbert-transform  derived 
signal-envelope  along  a  ray  is  depicted  as  an  example  of  input 
data  measured  on  a  phantom.  The  corresponding  attenuation 
profile  derived  from  the  data  is  presented  on  fig.  3  and  the 
signal,  restored  accordingly,  on  fig.  4. 

A  scan  of  a  pig’s  heart  (in-vivo)  is  shown  on  fig.  5, 
in  video  (envelope)  signal  representation.  The  attenuation 
map,  as  the  result  of  the  described  computation,  utilising  the 
same  data  as  those  used  for  the  image  fig. 5,  is  presented  on 
fig.  6  (higher  brightness  represents  a  higher  attenuation). 


Fig.  4  Restored  A-mode  signal 


Fig.  6  Computed  attenuation  map  of  the  object  on  fig.  5 


IV.  Conclusions 


References 


The  suggested  method  to  estimate  the  ultrasound 
attenuation  map  automatically  from  the  measured 
ultrasonographic  data  seems  to  provide  a  reasonable 
approximate  solution  to  the  given  problem,  thus  forming  a 
step  forward  on  the  way  to  hilly  automatic  attenuation 
compensation  in  ultrasonography.  The  model  of  ultrasound 
propagation  has  been  simplified  by  averaging  the  frequency 
dependence  and  also  by  neglecting  the  noise  effects,  then 
discretised  and  the  solution  sought  by  the  least-square-error 
method.  In  spite  of  the  mentioned  simplifications,  the  method 
has  proved  to  work  reasonably  well  for  available  sets  of  test 
data. 

Nevertheless,  the  independent  analysis  of  individual  ray 
A-mode  signals  does  not  take  into  account  the  lateral 
correlation  among  data  in  a  scan  (frame),  thus  producing 
some  artefacts  in  the  derived  maps:  the  relatively  large 
differences  between  neighbouring  A-mode  lines  are  in  some 
cases  rather  unrealistic.  Therefore,  the  two-dimensional 
analysis  of  the  problem,  taking  into  account  mutual 
dependencies  among  neighbouring  rays,  is  planned  as  the 
next  step  in  the  project. 
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